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Based on new developments in the understanding of supernovae (SNe) as gravitational- wave (GW) 
sources we estimate the GW background from all cosmic SNe. For a broad range of frequencies 
around 1 Hz, this background is crudely comparable to the GW background expected from standard 
inflationary models. While our estimate remains uncertain within several orders of magnitude, the 
SN GW background may become detectable by second-generation space-based interferometers such 
as the proposed Big-Bang Observatory (BBO). By the same token, the SN GWs may become a 
foreground for searches of the inflationary GWs, in particular for sub-Hz frequencies where the SN 
background is Gaussian and where the BBO will be most sensitive. SN simulations lasting far beyond 
the usual cutoff of about 1 s are needed for more robust predictions in the sub-Hz frequency band. 
An even larger GW background can arise from a hypothetical early population of massive stars, 
although their GW source strength as well as their abundance are currently poorly understood. 
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I. INTRODUCTION 

Core-collapse supernova (SN) explosions are among 
the most violent astrophysical phenomena. The total 
energy of about 3 x 10^'^ erg is primarily released in a 
neutrino burst lasting for a few seconds. About 1% goes 
into the actual explosion while only a fraction of about 
10~^ is emitted in visible light, yet a SN can outshine its 
host galaxy. In addition, SNe are expected to be strong 
gravitational-wave (GW) sources. About 1 SN per sec- 
ond takes place in the visible universe, but the neutrino 
burst has been observed only once from SN 1987A, the 
closest SN in modern history. GWs are even more dif- 
ficult to detect and have never been measured from any 
source. On the other hand, the next galactic SN proba- 
bly will be observed with high statistics in several large 
neutrino detectors as well as with the existing GW anten- 
nas. For the first time there is a realistic opportunity to 
observe such a cataclysmic event in all forms of radiation. 

However, galactic SNe are rare so that the next SN 
neutrinos to be observed could well be the diffuse neu- 
trino background from all cosmic SNe. The limit from 
Super-Kamiokande already touches the upper range 
of theoretical predictions HiSISl- The proposed gadolin- 
ium upgrade j5j or a corresponding future megatonne de- 
tector may well produce a positive detection. Therefore, 
it is natural to ask the analogous question whether the 
diffuse GW background from all past cosmic SNe will be 
important for current or future GW antennas that search 
for signals from individual astrophysical events or for the 
stochastic cosmic background that probably arises in the 
very early universe during the infiationary epoch. 

The existing literature focuses on GWs from the SN 
bounce signal in the kHz range that is relevant for bar 
detectors or ground-based laser interferometers such as 
LIGO, VIRGO, GEO and TAMA. However, recent stud- 
ies of SNe as GW sources indicate that a much stronger 



signal may be expected from the large-scale convective 
overturn that develops in the delayed explosion scenario 
during the epoch of shock-wave stagnation that may last 
for several 100 ms before the actual explosion 0, 
Moreover, the relevant frequencies reach below 1 Hz 
where they may be relevant for space-based detectors. 
While we find that a first-generation instrument such as 
LISA is not sufficient, the GW background from SNe may 
become relevant for a second-generation detector such as 
the proposed Big-Bang Observatory (BBO). 



While a detection of the GW background from SNe 
would be intriguing from the perspective of SN physics, 
particularly in conjunction with a future detection of 
the corresponding neutrino background, this possibility 
could be a significant problem for searches of the infia- 
tionary GW background that would probe a very early 
epoch of the universe. In fact, if ordinary astrophysi- 
cal foreground sources were seriously expected to mask 
the infiationary GWs in some range of frequencies, prob- 
ably an instrument like the BBO should be designed to 
have its optimal sensitivity in a different frequency range. 
Therefore, a reliable understanding of such foregrounds 
is an important task in view of the long-term perspective 
of observing the big bang in GWs. 



We begin our study in Sec. II with a general expres- 
sion for the GW background from all cosmic SNe. In 
Sec. HI we consider core-collapse SNe as sources, whereas 
in Sec. IV we turn to an early population of massive stars. 
In Sec. V we estimate the frequency range where the GW 
background is Gaussian before concluding in Sec. VI. We 
always use natural units with h = c = 1. 
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II. ESTIMATING THE GRAVITATIONAL 
WAVE BACKGROUND 

We assume all SNe are identical GW sources defined 
by the Fourier transform h{f) = dte~'^^'^f'^h{t) of 
the dimensionless strain amplitude h(t), which is pro- 
portional to deviations from spherical symmetry of the 
energy-momentum tensor. In general, h{t) is obtained 
from numerical simulations. However, the contribution 
of asymmetric neutrino emission is explicitly |^ 



h{t) 



2Gn 

~1T 



t-D 



dt' L,{t')q{t'), 



(1) 



where Gn is Newton's constant, D the distance to our 
"standard SN," and Li,{t) the neutrino luminosity. Fur- 
ther, q{t) < 1 is an asymmetry parameter defined as the 
angle-dependent neutrino luminosity folded over an an- 
gular function: for details see Eqs. (27)-(29) in Ref. 0. 

Equation shows that h{t) converges to a constant 
value hao for t }^ t,,, with t,, the neutrino emission time 
scale of a few seconds. Therefore, the neutrino burst 
causes h{t) to jump from zero to a non- vanishing constant 
value, an effect called the "burst with memory" 0, 
in 111 in • We now write 



-i^ / dtL,{t)q{t), (2) 



where in the second step we used Eq. Q and assumed 
/ ^ (27rtf;)~\ which is known as the zero-frequency 
limit Uliatlllll- This imphes 



27r 



(3) 



2.6 X 10"^^ {q) ^ 



lOkpc 
D 



3 X 10^3 erg 



where E^, = dt' L^{t') is the total emitted neu- 

trino energy and {q) = dt' L,^{t')q{t')/ E^ is the 

luminosity-weighted average neutrino anisotropy. Al- 
though the zero-frequency limit Eq. ||2Jl is valid for / <C 
{2TTte)~^, we will use it to evaluate the GW signal at 
/ ^ (27rtsim)~^, where tgim < is the maximal simula- 
tion time (less than a second) , and then we continuously 
extend to lower frequencies. Current simulations 

do not cover the strain spectrum below fractions of a 
Hertz. 

The energy density in GWs at frequency / per loga- 
rithmic frequency interval in units of the cosmic critical 
density pc = 3H^ / {SttG^^s) can be written as [T6| 



^gw(/) 



15GnPc 



dz 



ft 



(4) 



where Rsniz) is the SN rate per comoving volume, = 
{l + z)f, and h(f) is the strain spectrum of an individual 
SN at distance D. The cosmological model enters with 
\dt/dz\ = [(1 -I- z)H{z)]~^ and, for a flat geometry. 



H{z)^Ho [nMii + zf + riA] 



1/2 



(5) 



We will use the parameters Hm = 0.3, fl\ — 0.7, and 
Ho = ho 100 km s^^ Mpc"^ with ho = 0.72. 
The SN rate as seen from Earth is 



°°^^ i?SN(^) dV 



1 + z dz 



Aur^ ( z^ 
dzRs^{z) y , (6) 

{^ + z)H{z) 



where dV/dz is the fractional volume element, the cosmic 
expansion rate at redshift z is given by Eq. Q, and r{z) 
is the comoving coordinate, dr — (\ + z)dt. 



III. SUPERNOVAE 

The cosmic star-formation rate and, as a consequence, 
the core-collapse SN rate is reasonably well known at red- 
shifts z ^ 5 ,1^], but becomes rather uncertain at larger 
z. The SN rate per comoving volume is often expressed 
as 



Rsn{z) — Rs 



SN 



2'^-°'{l + zY 



for 2; < 1 , , 

for 1< z < 20 ■ 



If not otherwise stated, we will use jS = 2.7, a 



0, 



and the present-day rate Rg^ — 2 x 10~^ Mpc""^ yr~^. 
These values are consistent with the Super-K limits on 
the diffuse neutrino background The parameter ^ 
a < 2 is much less constrained than (3 and i?gj^ because 
a influences the rate only for z > I where the neutrinos 
are redshifted below the Super-K threshold This 
situation may change if low-threshold detectors such as 
the gadolinium upgrade of Super-K 5] become operative. 

For the strain spectrum h{f) of an individual event we 
first use the recent 3-dimensional asymmetric simulation 
shown in Fig. 9 of Ref. 0- We plot in Fig. the resulting 
GW background together with an estimate scaled down- 
ward by a factor of 100. The latter GW background with 
a total energy released in GWs of about 8 x 10""'^^ Mq 
should be considered more realistic . We also show sen- 
sitivities of LISA, BBO _18], LIGO (or EGO) correlated 
third generation (27(1, and the ultimate DECIGO 
which would be a quantum limited space-based interfer- 
ometer with mirror masses of 100 kg. For some frequen- 
cies, our lower estimate for the GW background from 
SNe can be comparable to the most optimistic GW back- 
ground from slow-roll inflation (horizontal lines) . 

Next, we use the three models sl5r, sllnrlSO, and 
pnslSO from recent two-dimensional (2D) core-collapse 
simulations based on a detailed implementation of spec- 
tral neutrino transport or following the long-time evo- 
lution of the newly born neutron star with good reso- 
lution, respectively 0. Models sl5r and sllnrlSO are 
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FIG. 1: Spectrum of GW background from Eq. |0J for a simu- 
lation from Ref. 7]. Solid lines: the spectrum is continuously 
extended with the low-/ tail according to Eq. (|5Jl. The up- 
per curve is obtained from Fig. 9 of Ref. i7|i whereas the 
lower is the same shifted downwards by a factor 100, which 
is considered a more realistic estimate. Horizontal lines: GW 
stochastic spectrum produced during slow-roll inflation and 
evaluated from Eq. (6) of Ref. We assume TjS = 0.3 

for the ratio of the tensorial and scalar contributions to the 
cosmic microwave background radiation (CMBR) anisotropy 
and consider two values ±10~^ of the running of the tenso- 
rial power-law index. Dotted lines: sensitivities of the indi- 
cated detectors. The BBO sensitivity is approximate and may 
change slightly in the final design. 



fully 2D, i.e. axially symmetric, stellar core collapse and 
SN simulations of a rotating 15 and a non-rotating 
11.2 Mq star with very strong core convection after shock 
formation. Convection takes place both in the neutrino- 
heated postshock layer and inside the nascent neutron 
star 0, |23| ■ The simulations were done with a full spec- 
tral treatment of neutrino transport and neutrino-matter 
interactions with gravitational rcdshift effects taken into 
account. An approximation to 2D transport was used 
which takes into account the dependence of the trans- 
port and neutrino-matter interactions on the laterally 
varying conditions in the star. It is based on the so- 
lution of the coupled set of radiation moments equa- 
tions and Boltzmann transport equation, assuming lo- 
cally radial neutrino fluxes. Moreover, order vjc veloc- 
ity terms in the moments equations due to the motion 
of the stellar fluid and the momentum transfer to the 
stellar plasma by neutrino pressure gradients were taken 
into account ( "ray- by-ray plus" approximation; for de- 
tails see Ref. |21|). Model pnslSO is a 2D hydrodynamic 
simulation of a non-rotating, convecting "naked" proto- 
neutron star, which was evolved through its neutrino- 
cooling phase from shortly after core bounce until about 
1.2 s later IH IH [23. The SN explosion was as- 
sumed to have taken off and therefore accretion and post- 
shock convection were not included in this model. Neu- 
trino transport was treated by flux-limited equilibrium 



diffusion, again assuming that the neutrino flux follows 
locally radial gradients according to the varying condi- 
tions in the 2D stellar background. Assuming neutrinos 
in equilibrium with the stellar medium also implied that 
effects of neutrino advection and "radiation compression" 
and the fluid acceleration by neutrino pressure gradients 
were accounted for. In all models the neutrino emission 
developed a time-dependent anisotropy because of the 
two-dimensionality of the evolving stellar background. 
The lateral variation leads to a non-vanishing value of 
The approximations in the description of 2D trans- 
port tend to yield an underestimation of the power in 
emission asymmetries with large angular scales and low 
frequency, but an overestimation of high-frequency and 
short- wavelength anisotropics in the neutrino emission. 
The reason for this behavior is that local fluctuations in 
the convective layers around the neutrinosphere trans- 
late by the mostly radial transport into small-scale fluc- 
tuations of the neutrino emission. In reality, however, 
multi-dimensional transport is much more non-local and 
the emitted radiation therefore contains averaged infor- 
mation from all parts of the emitting surface, thus aver- 
aging over short-wavelength variability of the hydrody- 
namic medium of the star (for a discussion in the context 
of neutrino transport in SN cores see Ref. [26]). 

We first focus specifically on model sl5r which gives 
more optimistic GW signals. Figures |21 and 13 show neu- 
trino luminosity L,j(t), anisotropy q{t), and GW strain 
h{t) and GW source spectra f\h{f)\, respectively. Fig- 
ure 131 shows that anisotropic neutrino emission domi- 
nates at low frequencies. Note that the simulation stops 
at about 250 ms after the bounce when only about 
0.5 X 10^^ erg has been emitted in neutrinos, i.e. about 1/6 
of the total, and when / dtL^{t)q(t) — 2.55 x 10^° erg. 
The average anisotropy over this time period is (q) ~ 
0.45%, but from Fig. |21we can see that \q{t)\ tends to 
increase. If q{t) were constant, according to Eq. J^l, 
h(t) would increase proportional to the emitted neutrino 
energy, whereas for a fluctuating q{t) we could expect 
that h{t) increases roughly as the square root of that en- 
ergy. Therefore, an estimate of the GW signal h{t) when 
3 X 10^'^ erg has been emitted in neutrinos, can be ob- 
tained introducing in the strain an enhancement factor 
between ~ \/6 and ~ 6. 

The resulting GW backgrounds are shown in Fig. 0] 
For model sl5r the blue band corresponds to the uncer- 
tainty of the enhancement factor discussed above, and 
the red band is a rough estimate of the uncertainty due 
to redshift evolution. For 100 Hz < / < 1 kHz the signal 
is dominated by the convective motion in the neutrino- 
heated postshock layer, whereas at lower frequencies it is 
mostly due to asymmetric neutrino emission. This sce- 
nario corresponds to a total energy release in GWs of 
~ 1.8 X 1O~^M0 for the most conservative enhancement 
factor. Note that, although the zero-frequency signal is 
lower, the GW energy release is higher than in the case 
of Fig. n because it is dominated by high frequencies. 

Equations and Q show that for our standard red- 
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FIG. 2: Neutrino luminosity Lv{t), anisotropy q{t), and GW 
strain h{t) times distance from anisotropic neutrino emission 
only, see Eq. Q, for model sl5r of Ref. 6], as functions of 
time after bounce. We also show the GW strain h(t) times 
distance from anisotropic neutrino emission, using the average 
anisotropy (g) — 0.45%. 
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FIG. 3: GW source spectra: Solid and dotted ragged lines are 
total and neutrino contribution, respectively, for model sl5r 
of Ref. Q at distance D = 10 kpc. The straight solid line is 
the low-/ tail Eq. (|3J with hoc from Eq. Q, using q{t) and 
L^{t) of Fig.H leading to f\h{f)\ ~ 2.21 x 10"^^ The dashed 
line is the schematic spectrum for PopIII stars of Eq. ((SJ. 

shift evolution (a = and [3 — 2.7) and i?,^ ~ 3 x 10^^ erg 
liberated in neutrinos, to completely mask the inflation 
background down to 10 niHz we would need an average 
anisotropy {q) ~ 6%, which is much higher than what is 
predicted by our most optimistic simulation model sl5r 
where (g) ~ 0.45%. We note that average quadrupole 
anisotropics {q) of the order of 1% are currently not 
excluded by astrophysical observations. In fact a 1% 
dipole anisotropy corresponds to observed kick velocities 
~ 400 kms""'^ for a neutron star of mass ~ 1.4M0, but 
we do not know whether and how quadrupole and dipole 
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FIG. 4: GW background for model sl5r with rotation (colored 
bands and solid lines) and model sllnrlSO without rotation 
(dashed line) from Ref. 0|. The source spectra have been 
continuously extended using the zero-/ tail, Eq. for / < 
1 Hz, except for the solid lines. The blue band and lower 
solid lines for model sl5r reflect the plausible range 6-36 of 
the enhancement factor to correct for the limited time of the 
simulation. For model sllnrlSO the lower enhancement factor 
of 6 was assumed. We always use q = in Eq. ((TJ, except for 
the red band and upper solid line which show the difference 
between a = and q = 2 (for enhancement factor 36) for 
model sl5r. 



anisotropics are correlated. 

For a < 2 in Eq. {7)), the integral in Eq. Q converges 
rather quickly. Cutting the integral at z = 5 instead of 
z — 20, for example, lowers the predictions by less than 
a factor 2. In contrast, the poorly constrained redshift 
evolution for z > 1 causes more significant uncertainties 
as illustrated in Fig. 0] 

To test the low-frequency dependence on the fluctua- 
tions of the neutrino anisotropy, we analyze GW signals 
for the proto-ncutron star model of Ref. 6] whose evolu- 
tion was followed for 1.2 s. In Fig.[5lwe compare the GW 
spectra for (i) the full simulation integrated over 1.2 s 
and (ii) the spectrum obtained from replacing the time- 
dependent anisotropy q(t) of the neutrino emission with 
its luminosity- weighted average over the 1.2 s duration 
of the simulation, (g) ^ 3 x 10^^. The huge difference of 
(q) for models sl5r and pnslSO arises from the rotation 
in model sl5r, that leads to a global, slowly growing de- 
formation and thus produces a systematic trend on q(t). 
In contrast, the merely convective fluctuations in model 
pnslSO on time scales of a few milliseconds cause a much 
lower time-averaged value of q. 

From Fig. Owe conclude the following. For / > 30 Hz, 
replacing q{t) by (q) makes no significant difference (the 
blue curve is hidden below the red one) because the sig- 
nal is dominated by convection and thus does not depend 
on q{t). In contrast, for 1 Hz < / < 10 Hz, fluctuations 
in q{t) increase the amplitude by a factor ~ 10. We 
also note the following technical point due to the prop- 
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FIG. 5: GW source spectra from mass motions plus aniso- 
tropic neutrino emission for the proto-neutron star model 
pnslSO of Ref. ^ at a distance of 10 kpc. The blue curve 
is for the full simulation with the fluctuating q{t) integrated 
over 1.2 s, whereas for the red curve q{t) was replaced with 
{(?> ~ 3 X 10-^ 
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FIG. 7: GW background resulting from the source spectrum 
for the proto-neutron star model pnslSO of Ref. |gl shown 
in Fig. 1^ without continuation to / < 1 Hz. We use a = 
in Eq. (|7|l and an enhancement factor of 3 (about 1/3 of 
the neutrino energy is emitted until simulation end). The 
solid blue line thus corresponds to the blue line in Fig. |3 
(fluctuating q{t)), whereas the solid red line is for q{t) replaced 
by {q) (red line in Fig. |5J . Other notation as in Fig. 0] 
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FIG. 6: Neutrino luminosity Lv{t), anisotropy q{t), and GW 
strain h(t) times distance from anisotropic neutrino emission 
only, see Eq. (Q, for model pnslSO of Ref. 0], as functions of 
time after bounce. We also show the GW strain h(t) times 
distance from anisotropic neutrino emission, using the average 
anisotropy {q) — 3 x 10~^. 



erties of Fourier transformations and the SN simulations 
used here: In the proto-neutron star model pnslSO the 
frequencies covered by the simulation do not extend to 
the regime where the zero- frequency limit applies, be- 
cause L^{t)q{t) oscillates with an almost constant ampli- 
tude over the entire duration of the simulation as seen 
from Fig. [SJ and thus has structure on time scales over 
which exp(i27r/t) still oscillates even in the lowest fre- 
quency bin. If neutrino emission would cease right after 
the end of the pnslSO simulation, fh{f) obtained with 



the fluctuating q{t), i.e. the blue curve in Fig. |S1 should 
approach f\h{f)\ 4 x 10^^"^ for / <C 1 Hz. In contrast, 
in model sl5r the product L^{t)q{t) starts out very small 
over more than half of the simulation time and tends to 
increase at later times, as seen from Fig. |5| Therefore, 
in the lowest frequency bin covered by the simulation of 
model sl5r, exp(i27r/t) does not fluctuate significantly 
over time scales on which Ly{t)q{t) has structure, and the 
zero-frequency limit overlaps with the lowest frequencies 
covered by the simulation. 

The GW background for the proto-neutron star models 
of Fig. [S] is shown in Fig. [T] Due to the redshift effect 
those GW spectra lie in part in the frequency range of the 
BBO space-based detectors. Here we are only interested 
in the relative amplitude between cases (i) and (ii), and 
not in their absolute amplitude, which is well below the 
sensitivity of space-based detectors. 

Finally, the "tail" of the neutrino luminosity during 
neutrino cooling (times later than 1-2 s after neutron 
star birth) is known from spherically symmetric neutron 
star cooling simulations [23, Is^, IMl and can be roughly 
approximated by L cx with n around unity (we use 
n = 1.1). This in principle could be used to extrapolate 
the signal to low frequencies. However, the evolution 
of the anisotropy parameter q{t) during the late phases 
of proto-neutron star cooling (> 1-2 s after bounce) has 
not been determined by models yet. Any assumption 
about its behavior after ~ 1 s implies major uncertain- 
ties. As discussed above, a constant anisotropy parame- 
ter or a fluctuating one can give very different predictions 
for the GW spectrum in different core-collapse models: 
For 1 Hz < / < lOHz, the GW source amplitude is sup- 
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pressed by roughly a factor 10 in model pnslSO if q{t) is 
replaced by a constant equal to its luminosity-weighted 
average (see Fig.El). In contrast, for model sl5r, the GW 
source amplitude is almost identical in these two cases in 
that same frequency range, recall Fig. |31 Similarly, if we 
evaluate the GW source spectrum by extrapolating the 
neutrino luminosity decline by L cx t~^-^ and replacing a 
fluctuating q{t) with a constant equal to its luminosity- 
weighted average, we would obtain an h{f) smaller by one 
order of magnitude in the frequency range of space-based 
detectors. Note that uncertainties in the GW source am- 
plitude of the order of 10 translate into uncertainties in 
figw of a factor ~ 100. 

In conclusion, temporal fluctuations of the neutrino 
anisotropy are extremely important for predicting the 
stochastic GW signal at frequencies / '--^ 0.1 Hz so that 
longer simulations are needed to make robust predictions. 



IV. POPULATION III STARS 

Core collapse events of the hypothetical population III 
(PopIII) generation of first stars [s^, 113 could be much 
more efficient emitters of GWs than today's SN popula- 
tions j34j . Lacking more detailed numerical models, we 
assume a schematic standard case of a mass of 300 Mq, 
with a shape of the GW spectrum at D = 10 kpc resem- 
bling that for ordinary SNe shown in Fig. |3I 
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16^2^2 



15G^ 



df\fh{f)\^ 



2 X lO-^M, 



(8) 



(9) 



for the total energy emitted in GWs. This very rough 
estimate is based on hydrodynamical mass motions in 
2D geometry If we use the zero-frequency limit of 
Eq. © and E^, ~ 10^^ erg emitted in neutrinos as given 
in Fig. 9 of Ref. [s^l for the model with rotation of their 
Table 1, we find that the average anisotropy would corre- 
spond to {q) ~ 3% in this case. While this is significantly 
larger than the values expected for ordinary SNe, it may 
not be unrealistic, given that PopIII stars tend to rotate 
more rapidly and explode more violently. 

Finally, we assume that the PopIII rate is concen- 
trated around a redshift zni — 15, which may explain 
reionization ll^. We thus write for the PopIII rate 
Rs^iz) = R6{z — ziii), where i? is a normalization con- 
stant. We further assume that all PopIII stars have a 
progenitor mass Mm = 3OOM0. The total number of 
PopIII events per comoving volume that occurred up to 
z = can then be written as 



dtRs-N{z) 



R 



(1 + zui)H{zin) 
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FIG. 8: GW background for our standard PopIII core collapse 



case and event rates 4.4 x 10 



< < 0.2 s- 



For 



the blue shaded band the source spectrum shown in Fig. |3] 
was cut off at / < 1 Hz and in this sense does not rely on 
the zero-frequency limit. The solid lines are obtained from 
continuing the source signal by the zero- frequency limit. 



where rij ~ 410 cm~^ is the CMBR photon density at 
redshift zero, rj ~ 6.3 x 10"^" is the number of baryons per 
CMBR photon, rriN is the nucleon mass, and /m is the 
fraction of all baryons going through PopIII stars. Elimi- 
nating R, it then follows that the rate observed at Earth, 
i?iii given in Eq. jnj, is directly proportional to /m. 



i?iii ~ 47rr2(zni) /iii%fy 



Mm 



(11) 



(10) 



For the parameters given above this yields Rm ~ 
0.2(/iii/10-3)s-i. 

The resulting GW background is shown in Fig. |H1 as 
a band delimited by two extreme assumptions about 
the total cosmic PopIII core-collapse rate observed from 
Earth, i.e. ii!iii ~ 0.2 s~^ (upper edge of band) and 
~ 4.4 X 10~* s~^ [l^ (lower edge). The zero-frequency 
limit would affect the GW spectrum only at frequencies 
below 0.05 Hz. We remark that certain models attempt- 
ing to explain the near infrared excess with PopIII stars 
require baryon fractions as high as /m ~ 0.05 cor- 
responding to rates i?iii ~ 10 s~^ and GW backgrounds 
a factor ~ 50 still higher than the upper band shown in 
Fig.|Sl Rates corresponding to Rm ~ 1 s~^ are discussed 
in Ref. [s^- The problem with such high-rate scenarios 
is to avoid overproduction of metals. 

Whereas neutrinos from PopIII stars are redshifted to 
energies probably too low to be detected 39], the GW 
background for 10 mHz ^ f 1 kHz could be dominated 
by these objects. 



V. GAUSSIANITY 

Finally, we estimate the duty cycle of the background 
at a given frequency /. It is defined as the product of 
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the rate of events and the time scale of interest, f~^. 

The duty cycle is roughly given by dividing the total 
cosmic SN rate, Eq. by the frequency. A probably 
more realistic estimate is obtained by folding the rate in- 
tegral over redshift with the redshift dependent contribu- 
tion to the GW signal given by the integrand in Eq. Q). 
For model sl5r of Ref. and the evolution discussed 
around Eq. Q these two estimates give approximately 
33(//Hz)-i and 4(//Hz)-i, respectively. 

The signal from ordinary SNe becomes Gaussian for 
/ 1 Hz, in agreement with a cosmic SN rate of about 
1 s^^. In contrast, the backgrounds in the 10^-10^ Hz 
range are not Gaussian, consistent with Ref. 0| where 
only the bounce signal was considered that lasts about 
10"'^ s. At these frequencies the duty factor is ^ 1%. 
The putative PopIII background becomes Gaussian only 
for / ^ ^iii, i-e. probably for / < 0.1 Hz. 

VI. CONCLUSIONS 

The GW background from cosmological core-collapse 
SNe at frequencies below 1 Hz is Gaussian. Its power 
is uncertain by several orders of magnitude, mostly due 
to uncertainties of SNe as GW sources. The most im- 
portant uncertain parameters are the asymmetry param- 
eter (q) of late-time neutrino emission and the param- 
eter a determining the star-formation rate at redshifts 
z > 1 [see Eq. Q]. The most optimistic current sim- 
ulations (for stellar core collapse with rotation) predict 
(q) ^ 0.45%, whereas Super-K limits on the diffuse neu- 
trino background constrain a in the range ^ a ^ 2. 
Using these parameter values and extrapolating the GW 
signal with the zero- frequency limit, Eq. ||2J), to the fre- 
quency band of space-based detectors, we showed that 
the GW background from SNe could be comparable to 
the maximum signal expected from standard inflationary 
models and thus could be detected by second-generation 
space-based detectors. However, to completely mask the 
inflationary background down to 10 mHz, we would need 
a very large asymmetry of (q) ^ 6%, that is not compati- 
ble with current simulations. By using the proto-neutron 
star model of Ref. , for which the simulation continues 
until 1.2 s, we observed that the GW spectra at frequen- 
cies of interest for future space-based detectors depends 
strongly on the fluctuations of q{t). Predictions could 
differ by a few orders of magnitude if q{t) continues to 
fluctuate for all or part of the neutrino cooling phase even 
after 1 s of neutron star cooling. 

Numerical simulations and stronger constraints or de- 
tections of the corresponding neutrino background will 
lead to an improved prediction of the GW background. 
However, accurate calculations of the GW emission from 
core collapse are hampered by our incomplete knowledge 
of important input physics like the equation of state at 
nuclear and supernuclear densities, of initial conditions 
like the structure of stellar cores at the onset of grav- 
itational collapse, and the amount of rotation and size 



of magnetic fields in the pre-collapse core. More reliable 
GW predictions also require a better understanding of 
the explosion mechanism and improvements in the nu- 
merical modeling of SNe. Quantitatively meaningful cal- 
culations will have to be done in 3D, will have to include 
a decent treatment of general relativistic effects, and will 
have to be done with full 3D neutrino transport in or- 
der to determine the neutrino emission anisotropy on all 
angular scales and at all frequencies. Three-dimensional 
effects are expected to reduce the emission anisotropy of 
neutrinos because of smaller convective structures |3 and 
because of the tendency of multi-dimensional transport 
to reduce the power on high spatial frequencies. In con- 
trast, significant rotation and large-scale magnetic fields 
might cause global emission anisotropics with possibly 
constant anisotropy parameter and high power at low 
frequencies. It is left to future generations of SN models 
to shed light on these important questions. 

PopIII stars could give a particularly strong contribu- 
tion, masking almost completely the inflationary back- 
ground, if we use the upper range of the very un- 
certain rate predictions |36l IstI l38j | in the literature and 
assume a total energy emitted in GWs of 2 x 10~^ Mq. 
For typical event rates and/or lower values of the to- 
tal GW energy emitted, the GW spectrum from PopIII 
stars could be comparable to or lower than the one ex- 
pected from inflation. In this case the zero-frequency 
limit would influence the GW spectrum only at frequen- 
cies below 0.05 Hz. In any event, the largest problem for 
a reliable prediction of the GW background from stellar 
core collapse is the almost complete lack of theoretical 
understanding of massive PopIII stars as GW sources 
and the extreme uncertainty of their abundance, in fact 
the uncertainty of their very existence. 

In summary, the uncertainties of the stochastic GW 
background from the core collapse of ordinary SNe and 
from PopIII stars is large, but it is intriguing that its 
power is in the ball park of the maximum GW back- 
ground expected from inflationary models. Of course, 
the inflationary GW background itself is also very uncer- 
tain and could be smaller than the benchmark shown in 
our figures. If stochastic GWs are detected with a future 
instrument like the BBC, it may be challenging to dis- 
entangle the inflationary GWs from those caused by SNe 
or the collapse of PopIII stars. 
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